Observing the spin Coulomb drag in spin-valve devices 
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Abstract 

The Coulomb interaction between electrons of opposite spin orientations in a metal or in a doped 
semiconductor results in a negative ofF-diagonal component of the electrical resistivity matrix - the 
so-called "spin-drag resistivity" . It is generally quite difficult to separate the spin-drag contribution 
from more conventional mechanisms of resistivity. In this paper I discuss two methods to accomplish 
this separation in a spin-valve device. 
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It is theoretically well established that the Coulomb interaction between electrons 

in a metal or in a doped semiconductor has a deeper effect on spin-polarized currents than on 
ordinary spin-unpolarized ones. The main reason for this is that the difference between the 
momenta of the up-spin and down-spin electrons is not conserved in a Coulomb scattering 
event: the transfer of momentum between electrons of opposite spin orientations therefore 
provides an intrinsic mechanism for the decay of a spin current, even in the absence of 

n 

electron-impurity scattering. This effect has been called "spin Coulomb drag" 1], or just 
spin drag for brevity. Mathematically, the spin-drag effect is best described in terms of the 
so-called spin transresistivity pn, which is defined as follows: Let j| and ji be the electrical 
currents associated with up- and down-spin electrons (we consider here for simplicity only 
currents in the x-direction and neglect spin-orbit effects), and let E^-, be the electro- 
chemical fields acting on the up- and down-spins respectively (The electro-chemical field E^- 
is defined as the gradient of the electro-chemical potential p^ divided by e. The electro- 
chemical potential itself is the sum of the true electric potential, which determines the 
position of the bottom of the conduction band, and the chemical potential, which determines 
the level of occupation of the band.) Then, for small departures from equilibrium one has 

Ea = ^ Paa'ja' , (1) 

n 

where the resistivity matrix po-o-' has the form p|| 

2 + 7r — 7r 
-7r 2 -|- 7r 

In the above equation 7 is the spin-drag coefficient, i.e., the intrinsic relaxation rate of 
the spin momentum — Pi, ^ is the ordinary momentum relaxation rate due to electron 
impurity interactions, m* and e are the band mass and the absolute value of the electron 
charge, and n is the total electronic density. Eq. (0) is valid under the assumption that the 
spin-flip scattering rate is negligible in comparison to 7 - a condition that should be well 
satisfied except at very low temperatures ^. We have also assumed, for simplicity, that the 
system is paramagnetic, i.e., n| = nj^ = |, so that = p^. Looking at Eq. Q we notice an 
important fact: pn is negative, because it takes a negative electro-chemical field to prevent 
an up-spin current from flowing when a down-spin current is present. On the other hand, 
the positivity of dissipation requires both eigenvalues of p to be positive - a condition that 
is obviously satisfied by Eq. ^ provided 7 is positive. 



P=^\ ~ ■ "I- (2) 
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FIG. 1: Spin- valve device for the measurement of the spin drag effect. A predominantly up- 
spin current is injected in a non-magnetic (NM) semiconductor via ferromagnetic (FM) electrodes 
between which a potential difference V is applied. The voltage probes are polarized oppositely to 
the injectors and therefore measure the down-spin electro-chemical potential. 

It is clear that an experimental determination of 7 would be of great interest since the 
value of this quantity is controlled by many-body correlations, which are intrinsic to the 
equilibrium state of the electron liquid. The main difficulty is that the spin transresistivity 
cancels out in the ordinary resistivity p = Bll^^^ so one has to devise an experiment that is 
somehow sensitive to the "spin resistivity" Pspin = ^^^2^^"^ ~ ~^ obvious way to 

proceed, first proposed in Ref. Q], is to measure the electro-chemical potential drop in the 
down-spin component when an up-spin-polarized current is driven into the semiconductor 
via highly spin-polarized ferromagnetic electrodes (spin injectors). In the limit that the level 
of spin-polarization p of the ferromagnetic electrodes is 100% "up" , and the spin diffusion 
length of the semiconductor is much larger than its geometrical length, the injected current 
is entirely in the up-spin component. Under these conditions, the electro-chemical potential 
"drop" for down-spins will be negative, if 7 is a finite positive quantity, and would vanish 
if 7 = 0. If, on the other hand, the polarization of the electrodes is less than 100% then 
the down-spin electro-chemical potential drop may remain positive, foiling our attempts to 
detect and measure 7. 

Thus, a very important question concerning this proposal is: how large should the polar- 
ization of the electrodes be so that one may observe a negative down-spin electro-chemical 
potential difference (as opposed to the trivially positive electric potential difference)? 

This paper is largely devoted to providing a sharp answer to this question. It will be 
shown (see Eqs. and (jlOj) below) that as p increases from to 100% there is a critical 
value of p, given by Pc = jr—, at which the electro-chemical potential drop for down-spins 



3 



switches from positive to negative. Measuring pc amounts therefore to a measurement of tt. 
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The experiment could be carried out in a three-layer spin-valve structure 
such as the one shown in Fig. 1. The two electrodes/spin injectors could be made out of a 
large-gf-factor II- VI semiconductor, e.g. Bei.Mnj,Zni_^,_ySe, where g ~ 100 12], which can 
be completely polarized by the application of a modest magnetic field. These electrodes are 
used to inject a spin-polarized current into a nonmagnetic (NM) lightly doped semiconductor 
(e.g., GaAs) and the total resistance across the electrodes is measured. The main physical 
assumptions underlying the proposed measurement are as follows: 

1. The spin drag effect is important only in the nonmagnetic semiconductor (GaAs), 
where the density of carriers is low. This is because it is theoretically well established 
that the spin drag increases in magnitude as the density of the electrons decreases [ll . 

2. The magnetic field, which is needed to spin-polarize the electrodes, has a negligible 
effect on the electronic states in the non-magnetic semiconductor, in which the gf-factor 
is small. 

3. The spin-resolved conductivities of the electrodes (t|| and scale in proportion to 
the corresponding electron densities, i.e., (J|| = ^^cr-^ and cr^ = ^-^a^ , where 

is the total conductivity of the homogeneous ferromagnet. (Of course itself may 
shghtly depend on p: this question will be discussed below.) 

The analysis is based on the equation for the electro-chemical potentials derived in Ref. 
In the one-dimensional geometry of Fig.l this takes the form 



where the 2x2 matrix Mg-o-' is, for our purposes, completely specified by its right eigenvectors, 
namely 

1 
1 

(the charge mode) with eigenvalue 0, and 



(4) 
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(the spin mode) with eigenvalue where L is the spin diffusion length. The solution of 
Eq. is straightforward. To make the best use of symmetry we assume that the semicon- 
ductor layer extends from x = —W/2 to x = W/2. The electro-chemical potentials are then 
odd functions of x [/iQ-(x) = — x)], and we can focus only on the region x < 0. In this 
region we write 
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X < 2 



f < X < 
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where J is the charge current, and a** are the conductivities of the electrodes and of 
the semiconductor, and L-^ and L"^ the spin diffusion lengths in the electrodes and in the 
semiconductor, respectively. Notice that the continuity of the charge current, J, is already 
built in Eq. (jH)). The three constants Co, Ci, and C2 are determined from the continuity 
of the two electro-chemical potentials and of the spin current — at x = —W/2. 

Their explicit forms are easily found to be 
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As mentioned above, the solution for x > is obtained by means of the symmetry relation 

/i^(x) = -/i^(-x). 

The behavior of the solution (expressed in units of eJW/a^) is shown in Fig. 2. Basically, 
we observe an accumulation of down-spin electrons and a corresponding depletion of up- 
spin electrons (i.e., /i| > at the left interface, where up-spin electrons are injected. 
The opposite happens at the left interface, where up-spin electrons are extracted. These 
spin accumulations effectively create a diffusion barrier, which increases the resistance and 
reduces the efficiency of spin-current injection. Under the assumption ^ W the electro- 
chemical fields, defined as the slopes of the elctrochemical potentials divided by e, are nearly 



FIG. 2: The behavior of the electrochemical potentials /if (dashed lines) and /i| (solid lines) 
calculated from Eq. © for 7T = and 7T = 1 and expressed in units of ^^p^ in the parallel- 
electrodes configuration. The semiconductor is in the region — ^ < < ^. The other parameters 
are p = 90%, ^ = 10, ^ = 2, and = 4. Notice the negative slope of H[{x) at x = when 
7r = 1: such a negative slope is an unmistakable signature of the spin drag. 

exactly uniform in the nonmagnetic region and their values are given by 



^ , , J JWp 

Notice that Ei is always smaller than E-^ and would tend to zero for p — 100% in 
the absence of the spin drag effect. This is because as the polarization of the electrodes 
increases, the down-spin component of the current must decrease: in the absence of spin 
drag this would imply that a gradient in chemical potential of down-spin electrons must be 
present to balance the electric field, resulting in Ej^ ^ 0. The spin drag upsets this balance. 
It is now necessary to have a finite, negative E^ in order to balance the momentum transfer 
from up- to down-spin electrons. The change in sign in Ei is an unmistakable signature of 
the spin Coulomb drag and occurs when the spin polarization of the electrodes exceeds the 
critical value 



'1 + 4a2 sinh^ (^) + 2a sinh (f ) ^ - 1 

Pc 



2a sinh (^) 

^ , (10) 

1 + 7^ 
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where a = (notice that > 1 for 7 = 0, that is, in the absence of spin drag. For 
7r = 1 with the parameters of Fig. 2 we have Pc ~ 0.85). Thus by measuring the value of 
p at which Ei changes sign one can determine 7r. 

The main drawback of such an experimental design (which is conceptually analogous to 
the design of the Coulomb drag measurement in bilayer systems j?!]) is the need to establish 
separate electrical contacts for the up- and down-spin electrons. This could be accomplished 
by the introduction of ferromagnetic voltage probes, polarized oppositely to the current 
leads. Unfortunately, such "probes" are technically difficult to implant and complicate the 
analysis of the experiment, for they disturb the equilibrium distribution of the spin in the 
semiconductor. For this reason I now describe what should be a simpler method to determine 
the quantity 1 + 7r. The idea of the measurement is simply to compare the total resistance 
R of the circuit at p = (i.e., for unpolarized electrodes) and p = 1 (i.e., 100% spin- 
polarized electrodes). No spin-polarized voltage probes are required. We assume that the 
homogeneous resistance of the electrodes and the external wires (denoted by Rc for brevity) 
is small compared to the resistance of the non-magnetic semiconductor. The polarization 
dependence of Rc presumably amounts to an even smaller correction. At p = the total 
resistance is thus essentially equal to the ordinary resistance of the semiconductor: the spin- 
drag effect is invisible here. At p = 1, on the other hand, the resistance depends very much 
on whether there is spin drag or not. If the spin drag were absent, then the resistance would 
be twice the ordinary resistance of the semiconductor, because only one of the two spin 
channels is open to conduction. In the presence of spin drag the flow of the up-spin current 
is hindered by collisions with down-spin electrons, which are stationary on the average: as a 
result, the resistance of the conductor becomes more than twice the ordinary resistance - in 
fact we will show that it is 2+'-jt times the ordinary resistance. Thus, by taking the difference 
R{p = 1) — R{p = 0) and dividing it by R{p = 0) we arrive at an experimental determination 
of jT. It should be noted that in taking the difference R{p = 1) — R{p = 0) the resistance 
Rc of the wires and the electrodes largely cancels out, except for its polarization-dependent 
component, which we feel justified in neglecting. Furthermore, this determination does not 
depend on the value of the spin diffusion length in the semiconductor, L^, provided the latter 
is much larger than the length of the semiconductor itself, W - a, condition that should not 
be too difficult to satisfy in practice. Likewise, the value of the spin-diffusion length in the 
electrodes, L-^, is essentially irrelevant as long as the potential drop is measured between 
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FIG. 3: Behavior of the parahel-electrodes magneto-resistance Rparaiieiijp) — Rc{p) (in units of the 
ordinary resistance of the semiconductor, Rg) as a function oi p for 7r = 0, 1, and 2. Notice the 
sharp enhancement caused by the spin-drag resistivity for values of p close to 100%. As explained 
in the text, this can be used to determine 7T. 

points that are much farther than a distance from the FM-NM interfaces. 

Far from the FM-NM interfaces (more precisely, at a distance much larger than L^) 
the electro-chemical potentials of the two spin orientations tend to common values ^- = 



'JWiCo:zll + £^ for a; ^ and fi+ = ^ ejx ^ _^ rj.^^ difference 



between these two asymptotic behaviors is e times the voltage drop due to the presence of 
the semiconductor layer. Hence, the resistance of our device (per unit cross-sectional area) 
is given by 



where Rcip) is the combined resistance of the electrodes and the external wires, and W = ^ 
is the ordinary resistance of the semiconductor. The last term on the right hand side of this 
equation arises from the spin accumulations at the interfaces between the electrodes and 
the semiconductor. Fig. 3 shows the behavior of the key quantity Rparaiieiip) ~ Rc{p) as a 
function of p. It increases from i?s at p = to i?*(2 + 7r) at p = 1. Interestingly most 
of the change occurs in the region of p close to 1. This can be exploited to reduce the 
undesired effect of the p-dependence of Re- Namely, rather than considering the change in 
resistance from p = to p = 1, it may be sufficient to consider the change from say p = 0.5 
to p = 1.0 with correspondingly less variation in R^. Notice that heoretical calculations of 
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7 as a function of temperature and electronic density can be found in Refs. The 
temperature dependence of 7 is particularly interesting as it exhibits a characteristic broad 
maximum at about the degeneracy temperature of the carriers in the semiconductor. 

For completeness, let us now see what happens in the antiparallel-electrodes configuration. 
In this case, the electro-chemical potentials obey the symmetry relation fia{x) = —f^a{~x) 
and it is easy to see that the new solution is now obtained from the parallel case solution 
simply by interchanging the quantities sinh (^) and cosh (^). More precisely, the solution 
for X < takes the form 




where the constants Cq — C'2 are given by 
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(13) 
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The solution for x > is obtained by means of the symmetry relation (Xa{x) = — /x_ 
A representative plot of /i| and /i^ is shown in Fig. 4. 

We can calculate the resistance of the antiparallel-electrodes configuration in precisely 
the same way as in the parallel-electrodes case. The result is 

' W " 



2Wp'^ 

Rantiparallel = Rc{p) + + „, r siuh 



2L' 



(15) 



and the quantity Rantiparaiieiip) — Rc{p) is plotted vs p in Fig. 4. The resistance of this 
configuration is of course much larger than that of the parallel configuration (this is the 
well known GMR effect) and it is easy to see that in the limit p — > 1 it tends to to i?''(l + 
7r) (^)^. Notice that the Coulomb enhancement in this configuration is very sharply 
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FIG. 4: The behavior of the electrochemical potentials /i| (dashed lines) and /i| (solid lines) 
calculated from Eq. H12|) for = and 7r = 1 and expressed in units of ^^p^ in the antiparallel- 
electrodes configuration. The parameters are the same as in the caption of Fig. 2, namely p = 90%, 
^ = 10,^ = 2, and ^ = 4. 
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FIG. 5: Behavior of the antiparallel-electrodes magneto-resistance Rparaiiei{p) ~ Rc{p) (in units of 
the ordinary resistance of the semiconductor, Rg) as a function of p for 7r = 0, 1, and 2. The 
enhancement caused by the spin-drag resistivity for values of p close to 100% is now amplified by 
a factor (^r-)^ (= 16 in the present case). As explained in the text, this can be used to mesaure 
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once 7r is known. 
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confined to the region of j9 ~ 1. Tfie results of the above calculation can be used to 
determine once •jr has been determined from the measurement of the resistance in the 
parallel-electrodes configuration. 

In summary I have theoretically analyzed in this paper two methods to measure the 
spin drag coefficient of a non-magnetic semiconductor sandwiched between highly spin- 
polarized ferromagnetic electrodes. The first method builds upon the gedanken experiment 
proposed in Ref. |l| showing that an unambiguous qualitative signature of the spin drag 
effect occurs when the spin polarization of the ferromagnetic electrodes exceed the critical 
value Pc — Y+i^- second, more quantitative method one simply measures the extra 

resistance introduced by the relative motion of the up-spin and down-spin electrons in the 
semiconductor region of a basic spin-valve device. It is hoped that these discussions will 
encourage further experimental work aimed at the observation of the spin Coulomb drag. 
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